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Using numerical techniques, the diagonal and off-diagonal superconducting one-electron Green's 
functions are calculated for a two-dimensional (2D) t-J model on a periodic 32-site cluster at low 
doping. From these Green's functions, the momentum and frequency dependence of the pairing gap 
A(k,w) are extracted. It has d x 2_ y 2 symmetry and exhibits ^-dependent structure which depend 
upon J/t. We find that the pairing gap persists down to small J/t values. The frequency- and 
momentum-dependent renormalized energy and renormalization factor are also calculated. 

PACS numbers: 75.10.-b 71.27.+a 75.50.Ee 75.40.Mg 



For the last decade, the search for superconduc- 
tivity in models of strongly-correlated fermions, has 
been triggered by accumulating experimental evidence 
in favor of an unconventional (i.e. not based on the 
usual phonon-mediated interaction) mechanism in high- 
T c two-dimensional (2D) or ladder || superconduct- 
ing cuprates. Soon after the discovery of the layered 
cuprates, it was proposed that the Hubbard model || 
and its strong coupling limit, the t-J model captured 
the generic features of these materials. Early RPA ||, |(| 
and Gutzwiller variational calculations |Q supported the 
notion that the doped state would have d x 2_ y 2 pairing. 
Nevertheless, despite the conceptual simplicity of these 
Hamiltonians and the results of these approximate calcu- 
lations, the nature of the basic mechanism responsible for 
pairing as well as the actual physical properties of these 
models remained controvertial. In recent years, numeri- 
cal calculations have provided insight into the second of 
these questions regarding the actual properties of these 
models. Specifically, Monte Carlo calculations showed 
that the ground states of the undoped models had long- 
range antiferromagnetic order ^ Exact diagonaliza- 
tion studies showed that holes doped into a t-J cluster 
can form d x 2_ y 2 pairs ]l0| , p"T| . In addition, numerical 
calculations support the view that the doped system can 
have low-lying stripe domain wall states |12| which can 
be stabilized by lattice anisotropies Jl3|, |14[ . While the 
static striped phase competes with superconductivity, it 
has been argued that the addition of a next-near-neighbor 
hopping fllq| or the use of periodic boundary conditions 
fl6f tips the balance in favor of the d x 2_ y 2 pairing phase. 
The point is that a variety of numerical calculations pro- 
vide evidence that indeed these models do exhibit the 
basic properties seen in the cuprates. However, in spite 
of this progress, the question regarding the nature of 
the basic mechanism responsible for the pairing remains 



open. Within the traditional BCS framework, one would 
look for a reflection of the pairing interaction by examin- 
ing the momentum and frequency dependence of the gap 
A(k, oj). Because of the relatively short coherence length 
and the relatively high energy scales J and t, one can hope 
to learn about both the momentum and frequency depen- 
dence of A(k, lo) from an exact diagonalization study of 
a cluster. Here we present a numerical study aimed at 
doing this for the t-J model and obtain the first results 
for the w-dependence of the gap. 
The 2D t-J Hamiltonian reads, 

H = Jy(Si • S - -riirij) 

+ Y,t ij {c{ (T c j , a + h.c.) (1) 

id," 

where the exchange integrals Jy and hopping terms are 
(for simplicity) restricted to nearest neighbor (NN) sites 
(called hereafter J and t) and c\ are projected fermion 

creation operators defined as c\ a (l — n^-o-). Hereafter, 
unless specified otherwise, t sets the energy scale. The 
conclusions drawn in this Letter are supported by Exact 
Diagonalisation (ED) studies performed on the square 
cluster of N=32 sites depicted in Fig. |l](a) slightly doped 
with up to 2 holes [jl7j . This finite-size system is partic- 
ularly appealing since it exhibits the full local symme- 
tries of the underlying square lattice as well as all the 
most symmetric k-points in reciprocal space as seen in 
Fig. |l|(b). Note that, although the hole doping is quite 
small, the hole occupation shown in Fig. [l](b) for a few 
k-points is rather consistent with a "large" Fermi surface. 

A superconducting ground state (GS) is characterized 
by Gorkov's off-diagonal one-electron time-ordered Green 
function 

F k (t) = - (Tc_ k ,_ CT (V2)c klCT (-t/2)) . ( 2 ) 
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FIG. 1: (a) Periodic square 32-site cluster; (b) Reciprocal 
space of the 32-site cluster. The area of the circles (drawn 
only for a few k-points especially relevant to our analysis) 
correspond to the hole occupancies (ck^c^ CT ) in the 2-hole 
GS at J = 0.3 (whose exact values are also shown on the 
plot). 



Close to half-filling, in a finite system, can be com- 
puted from, 



F k (z) = (N- 2|5_ k , 



1 



z - H + E 



-c Krj \N 



N-l 



(3) 



defined for all complex z (with Imz ^ 0). Here the num- 
ber of particles in the initial |JV) (half-filling) and final 
|7V — 2) (two- hole doped) GS differ by two, reflecting 
charge fluctuations in a SC state. For convenience, the 
energy reference E^-i is defined as the average between 
the GS energies of \N) and \N - 2). Note that these 
states are both spin singlets. In addition, they exhibit 
s-wave and d x 2_ y 2 orbital symmetries (both of even par- 
ity) respectively. Owing to these special features, it is 
straightforward to show that the analytical continuation 
of F\c(uj) to the real frequency axis is even in frequency 
and can be expressed as, 



F(k,w) = F^lu + ie) + F k (- 



ie) 



(4) 



where e is a small imaginary part. The superconduct- 
ing frequency-dependent gap function A(k, lu) is directly 
proportional to F(k, lu) and a simple analysis shows that 
both functions have d x 2_ y 2 orbital symmetry. In partic- 
ular, they identically vanish on the diagonals k x = ±k y . 

Pairing between holes should also be reflected in the 
structure of the (time-ordered) diagonal Green function 
G(k, lu). Here, it is convenient to define the finite size 
G(k, lu) as the sum of the following electron- (i.e. occu- 
pied states for lu < 0) and hole-like (i.e. empty states for 
lu > 0) parts, 
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FIG. 2: Diagonal and off-diagonal one-electron (time- 
ordered) Green functions versus frequency (in units of t) ob- 
tained on the cluster of Fig. |l| and an average hole doping of 
3%; (a) -ImG(k,u); (b) real and imaginary parts of F{\a, u) 
for e = 0.05; (c,d) comparison between the data of ReF and 
— ImF obtained with e = 0.05 and with e — 0.005. 



so that both Green functions F and G have the same set 
of energy poles. Here, the well-known continued-fraction 
method used to compute diagonal correlation functions 
such as G(k, lu) is extended to deal with off-diagonal ones 
such as F(k, lu) |9). Data are shown in Fig. ||(a- 
d). The spectral density Im G(k, lu) exhibits sharp quasi- 
particle-like peaks both above and below the Fermi level 
(lu = 0). Contrary to G(k, lu), F(k, u>) has a significant 
amplitude only at low energy, typically for \lu\ < 4 J as 
seen in Fig. |^(b) , reflecting the energy scale of the pairing 
interaction. Note that, due to the discreteness of the 
low energy spectrum of the cluster, a finite value of e 
is necessary to wipe out the irrelevant fast oscillations 
of the F (or G) Green functions (see e.g. Fig. ||(c-d)). 
However, the static limit lu — ► (together with e — > 0) is 
perfectly controled and has a physical meaning. 

The frequency-dependent gap function A(k, lu) can be 
extracted from the knowledge of the Green functions by 
assuming generic forms of a SC GS p0| at low energies, 



G(k,w) = z(k,w) 
F{k,Lu) = z(k,w) 



lu + e(k, lu) 



lu 2 - {e(k,Lu) 2 + A(k,w) 2 ) + te 
A(k,Lu) 

lu 2 - (e(k,Lu) 2 + A(k,w) 2 ) + It 



(6) 



where z(k, lu) and e(k, lu) are the inverse renormaliza- 
tion parameter and the renormalized energy (contain- 
ing the energy shift) respectively. In the BCS limit, 
the renormalized energy e(k, lu) = £k and the gap func- 
tion A(k, lu) = Ak do not depend on frequency (and 
z(k,Lu) = 1) so that F(k, to) exhibits only two poles at 
±(e k + A 2 .) 1 / 2 . However, it is necessary to assume an 
explicitely frequency dependence of the gap function in 
order to reproduce the secondary peaks seen in Fig. 0(d). 
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FIG. 3: (a) Real and imaginary components of the frequency- 
dependent gap function (in units of t) at low energies, (b) 
Dynamical spin structure factor in the 2 hole-doped 32-site 
cluster. The positions in energy of the first pole (magnon) and 
of the mean value of the spectrum are indicated by arrows. 



Using the parity of e(k, ui) and A(k, ui) in the vicinity of 
the real frequency axis one gets, 



A(k,u) 



2wF(k,w) 



G(k,uj)~G(k,-Lu) 



(7) 



Thus, from a numerical calculation of G(k, ui) and 
F(k, ui), one can obtain A(k, u>). 

Results for A(k, ui) with k = (w, 0) at the M point, 
are shown in Fig. 3(a). Despite the limited resolution 
in frequency of the ED data, A(k, ui) exhibits dynamic 
structure on a scale of energies several times J. We believe 
that finite-size effects have increased the onset frequency 
wa and that the corresponding time 2tt/u>\ should be 
viewed as a lower bound of the characteristic time scale 
of the pairing interaction. For comparison, the dynamic 
spin structure factor Imx(<\, ui) for the 2-hole doped 32- 
site cluster with q = (w, tt) is shown in Fig. 3(b). 

The quasi-particle gap at the gap edge, where e(k, 0) = 
0, is determined from the self-consistent solution of 



ReA(k,u = A k ) = A k 



(8) 



As shown in Fig. 3(a), A(k, ui) is real and essential con- 
stant over an energy region larger than the gap. Fur- 
thermore, as we will see, e(k, 0) for k = (tt, 0) is small 
so that Ak is to a good approximation given by A(k, 0). 
For k = (7r,0), the zero frequency limit of the gap func- 
tion A(k, ui — > 0) (which is purely real) computed by 
numerically taking the u> and e — > limits of eq. (]?]) |2l|] 
is plotted versus J in Fig. 4(a). In Fig. 4(b) the static 
u> = value of the renormalization energy e(k, 0), the 
renormalization factor z(k, ui — 0), and the equal-time 
pair amplitude (c^ k -o-cj^) for k = (tt, 0) versus J are 
also shown [Q. The values of A(fc,0) and z(k,0) are in 
good agreement with the results Ohta, et. al jl8| found 
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FIG. 4: (a) Zero-frequency gap and typical frequency da (di- 
vided by 4) of the gap function (in units of t) plotted vs J. 
For comparison, the magnon energy and the mean-value uii 
of the spin structure factor (see Fig. 3(b) is also shown, (b) 
Equal-time pair amplitude (cl k „ .cj c J) calculated at momen- 
tum k = (tt, 0) vs J/t. The renormalized energy e(k, 0) (in 
units of t) and the renormalization parameter z(k,ui = 0) are 
also shown as filled squares and triangles, respectively. 



by fitting their numerical t-J 4x4 and V18 x yl8 clus- 
ter calculations to a BCS quasi-particle form in which 
the frequency dependence of A and z were neglected. 
They found that the magnitude of the gap at the anti- 
mode varies as 0.25J to 0.3J at low doping and as seen in 
Fig. 4(a) we find A(/c, 0) ~ 0.3J. It is interesting to note 
that the data for A(fc,0) in Fig. 4(a) do not show any 
lower bound of J in contrast to the pair binding energy 
]To| , [Tl| . The spin structure factor is plotted in Fig. 3(b) 
showing a low energy peak (magnon) and a higher energy 
background that could be connected to the structures in 
the gap function at an energy ~ J. For comparison, 
the magnon energy and the mean- value ]23|] of the spin 
spectral weight are plotted in Fig. 4(a) together with the 
zero-frequency gap and the onset frequency wa in the 
gap function. While the similarities in the J-dependence 
of these quantities suggest that the dynamics underlying 
the pairing mechanism is related to the spin fluctuations, 
further work on other clusters, such as 2-leg ladders, are 
needed to sort out the influence of finite-size effects. 

We conclude this investigation with a further discus- 
sion of the parameters which enter eq. (^). Note that 
since these parameters depend upon both k and ui, the 
static ui — results which we will discuss are near the 
energy shell only for k values near the fermi surface and 
for A k small compared to the characteristic frequency 
variations set by J. In fact, even trying to select k-points 
near the "Fermi surface" defined by e(k, 0) = requires 
some care since the Fermi surface (FS) cannot be exactly 
defined on a finite cluster, especially for a small num- 
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TABLE I: Static limit (cj — 0) of the gap function, the renor- 
malized energy and the renormalization parameter at a few 
k-points away from the nodes and for J ft = 0.3. 
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k-points 


A(k,0) 


c(k,0) 


z(k,0) 


(tt/2,0) 

(tt,0) 
(7r,7r/2) 


0.064958 
0.095641 
0.186100 


0.20260 
0.11867 
0.08895 


0.21221 
0.13542 
0.12855 
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ber of holes. However, the smooth variation of the hole 
occupancy in the BZ [see Fig. 1(b)] suggests that this 
system should still pick up features of a slightly doped 
AF with a large FS. With this in mind, the static limit 
of the renormalized energy e(k, 0) and the renormaliza- 
tion parameter z(k, 0) were obtained from the calculated 
E,u -> limits of G(k,0), F(k,0), and A(k,0) by using 
eq. (||) for u> = 0. These are tabulated for values of k 
both near the nominal FS and away from it in Table 1. 
The renormalized energy e(k, 0) has roughly the energy 



scak uf Llit; 



quasi-pai Ucle-like eneigy pules at Ihe avail- 
able' iiiuiiieiiLa uf Lhe BZ. The small value uf t(k, 0) given 

by our analysis for k = (7r, 0) show that the excitation 
that we are probing is relatively close to the FS (defined 
here as having e(k, 0) = 0. Note also that the weight 
z(k, 0) of these excitations is small (« J/2) compared 
to the large incoherent background [see Fig. 2(a)]. Its 
J-dependence plotted in Fig. 4(b) is consistent with the 
J a power-law behavior found previously in earlier calcu- 
lations of a single hole propagating in an AF background 
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